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An analogy of the Fokker-Planck equation (FPE) with the Schrbdinger equation allows us to use quantum me- 
chanics technique to find the analytical solution of the FPE in a number of cases. However, previous studies have 
been limited to the Schrodinger potential with a discrete eigenvalue spectrum. Here, we will show how this ap- 
proach can be also applied to a mixed eigenvalue spectrum with bounded and free states. We solve the FPE with 
boundaries located at x = ±L/2 and take the limit i — ► oo, considering the examples with constant Schrodinger 
potential and with Poschl-Teller potential. An oversimplified approach was proposed earlier by M.T. Araujo and 
E. Drigo Filho. A detailed investigation of the two examples shows that the correct solution, obtained in this 
paper, is consistent with the expected Fokker-Planck dynamics. 
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1. Introduction 

The one-dimensional Fokker-Planck equation (FPE) for the probability density p[x, t), depending on 
variable x and time t, assumes the generic form (Jf" 



dpix.t) d . id 2 \D{x,t) 



p[x, t) 



(1.1) 



Here, the drift coefficient or force f(x, t) and the diffusion coefficient D{x, t) depend on x and t in general. 
The Fokker-Planck equation is related to the Smoluchowski equation. Starting with pioneering works by 
Marian Smoluchowski |1, 2], these equations have been historically used to describe the Brownian-like 
motion of particles. The Smoluchowski equation describes the high-friction limit, whereas the Fokker- 
Planck equation refers to the general case. 

The FPE provides a very useful tool for modelling a wide variety of stochastic phenomena arising 
in physics, chemistry, biology, finance, traffic flow, etc. |3-6]. Given the importance of the Fokker-Planck 
equation, different analytical and numerical methods have been proposed for its solution. As it is well 
known, the stationary solution of FPE can be given in a closed form if the condition of a detailed balance 
holds. The study of the time-dependent solution is a much more complicated problem. The FPE <1.1> with 
a general time-dependence and a special x-dependence of the drift and diffusion coefficients has been 
studied analytically in |7|] using Lie algebra. This method is applicable when the Fokker-Planck equation 
has a definite algebraic structure, which makes it possible to employ the Lie algebra and the Wei-Norman 
theorem. Generally, there are only a few exactly solvable cases. A simple example is a system with con- 
stant diffusion coefficient and harmonic interaction of the form f{x) - -dV{x)/dx with harmonic poten- 
tial V[x) ~ x 2 . The case with double-well potential is already quite non-trivial and requires a numerical 
approach |8]. 

The known relation between the Fokker-Planck equation and the Schrodinger equation can also be 
used. This approach allows us to apply the well known methods of quantum mechanics. In particular, 
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analytical solutions can be found in the cases, where the eigenvalues and eigenfunctions for the consid- 
ered Schrodinger potential are known. For a general Schrodinger potential, numerical treatments used 
in quantum mechanics, such as the Crank-Nicolson time propagation with implicit Numerov's method 
for second order derivatives 1 9], are very useful. To apply it to Schrodinger-type equation, we just need to 
replace the real time step A? by an imaginary time step Af — «■ -iAf. In quantum mechanics, this is called 
imaginary time propagation and is used for calculation of both ground states and excited states. The an- 
alytical studies of mapping the FPE to Schrodinger equation have been so far restricted to a treatment 
of discrete eigenstates. An attempt has been made in 1 10] to extend this approach to the potentials with 
a mixed (discrete and continuous) eigenvalue spectrum. However, we have found a basic error in this 
treatment, indicated explicitly in the end of section l4~3l 

The aim of our work is to show how the problem with mixed eigenvalue spectrum can be treated 
correctly. We will show this in two examples: one with constant Schrodinger potential and another with 
Poschl-Teller potential. The same example has been incorrectly treated in 1 10] . To avoid any confusion 
one has to note that the Poschl-Teller potential is referred to as Rosen-Morse potential in 11011 . 



2. Solution of FPE with constant diffusion coefficient 

We start our consideration with the one-dimensional Fokker-Planck equation dl.H in the following 
formulation 

dp{x,t) d r „ . Dd 2 p{x,t) 

for the probability density distribution p{x, t), depending on the variable x and time t. Here, f{x) is the 
nonlinear force and D is the diffusion coefficient, which is now assumed to be constant. We consider 
natural boundary conditions 

dp{x, t) 

lim p{x,t)= lim p ' — = (2.2) 

x—> +oo x— ±oo OX 

and take the most frequently used initial condition 

p(x,t = 0) =6{x-Xo) (2.3) 



in the form of the 5-function. This FPE {2.1) can be transformed into an equation of Schrodinger type 
(see section \2~2) . Unfortunately, the well known relation [see equation d2.25H , derived for the discrete 
eigenvalue spectrum, cannot be applied if this equation has a continuous or mixed eigenvalue spectrum. 
To overcome this problem, we follow a properly corrected treatment of 1 10]. Namely, we solve the FPE 
with boundaries located at x = +LI2 and then take the limit L^oo (see section lPV This approach is used 
in quantum mechanics to describe unbounded states. To keep a closer touch with quantum mechanics, 
here we will use the boundary conditions p{x — +LI2, t) = 0, further referred to as absorbing boundaries. 

2.1. The stationary solution 

The stationary solution p s t{x) is the long-time limit of p[x, t) at t — > oo, which follows from the equa- 
tion 

d , , Dd 2 p st (x) 

0=-[/ W p stW ]--^. (2.4) 

The force f{x) can be expressed in terms of the potential V{x) via f{x) = -dV{x)ldx. It yields 

(2.5) 



o- d 

dx 



dV{x) Ddp st {x) 

Pst(x) + 



dx 2 dx 

Due to the natural boundary conditions, we have zero flux 



M(x) = — PstW-- P ! =C with C=0. (2.6) 

dx 2 dx 
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Thus, we have 



dpst(x) = 2 dV{x) 
dx D dx 

dp s t(x) 2 , 

dV(x) , 



PstW , 



D 



which yields the stationary solution 
where 



Pst(x) 

p st {x) = Jf- 1 Y{x), 



Y{x) = exp 



(2.7) 
(2.8) 

(2.9) 
(2.10) 



has the meaning of an unnormalized stationary solution only in case of natural boundaries and JV is the 
normalization constant 

+oo 

2 



{ 



dxexp 



D 



V(x) 



(2.11) 



This function Y{x) is further used to construct a time-dependent solution. 



2.2. The time-dependent solution with discrete eigenvalues 

Here, we derive a time-dependent solution, starting with the transformation p{x, t) -* q(x, t) defined 

2 V{x) 



by 



p(x, t) = Y ll2 {x) q(x, t) = exp 



D 2 



q{x, t) . 



(2.12) 



This transformation removes the first derivative in the original Fokker-Planck equation and generates 
the equation of Schrbdinger type for the function q[x, t), i. e., 



dq{x,t) Dd 2 q{x,t) 
= -VsM4(x,f) + - 



where 



dt 



,ld 2 l/(x) 2 



dx 2 



1 dV{x) 



2 dx 



(2.13) 



(2.14) 



is the so-called Schrbdinger potential. In the case of discrete eigenvalues, we apply the superposition 
ansatz 

oo 

q(x,t) = Y,a„{t)ilf n (x). (2.15) 

n=0 



After inserting H2.15) into (2.13) , we get the eigenvalue problem 

D d 2 w n {x) 

~2dt 2 ^W?nW = -A„? B W (2.16) 

for eigenfunctions f n (x) and eigenvalues A„^0 with time-dependent coefficients a n {t) given by 

a n {t) = a„(0)exp(-A„ t) . (2.17) 
According to this, equation (2.15) can be written as 



q{x,t)= ^ a n {0)e *">„(*). 

n=0 



(2.18) 



The eigenfunctions i[/ n {x) are orthonormal, i. e., 

+oo 

J ■y n (x)yj m (x)dx = 8 nt 



(2.19) 
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and satisfy the closure condition (completeness relation) 



£lMxV,,(x) = $(*-*'). (2.20) 
n=0 



Equation <2.16> can be written as a Schrodinger-type eigenvalue equation with Hermitian Hamilton op 
erator 36 : 

D d 2 
2 dx 

The coefficients a n (0) in <2.18) are calculated using the initial condition 



jey n {x) = X n f n {x) with J£> = - — — + V${x) . (2.21) 



p(x, t = 0) = Y ll2 (x)q(x, t = 0) = 8{x- x ) . (2.22) 
According to 1 12. 18> . this relation can be written as 

oo 

y- 1/2 M5(x-x )= L a m {ff)fmix). (2.23) 

m=0 

In the following, we multiply both sides of this equation by y/ n {x) and integrate over x from -oo to +oo. 
Taking into account I2.19K it yields the so far unknown coefficients 

fl n (0)= Y- V2 {x )y/ n {x ). (2.24) 

The final result of this calculation reads 



/ Y(x) 00 

P {x ' t] = \ ~v~i — 7 E e V»CJ£b)V»W ■ (2.25) 
V * ( x oi „=0 

Note that this method can also be used for other boundary conditions. The solution in the general form 
of H2.251 is well known from older studies, e. g., 1 11] and can be found in many textbooks, e. g., 

2.3. The time-dependent solution with mixed eigenvalue spectrum 

Consider now the problem with two absorbing boundaries located at x — +L/2 instead of the natural 
boundary conditions. In this case, we have a discrete eigenvalue spectrum, and equation 12.25) can be 
used (with summation over exclusively those eigenfunctions which satisfy the boundary conditions in a 
box of length L) to calculate the probability distribution piix, t), i. e., 



/ Y(x) 00 

p L {x, t) = J —— £ e~ A - L t yr n ,axo)Vfn,ax) , (2.26) 

where X n> i are eigenvalues and ^(jr) are the corresponding eigenfunctions, which fulfill the boundary 
conditions. Let us split this infinite sum into two parts: for X n> i < A C on and X n> i 3= A con , where A con is the 
smallest continuum eigenvalue in the case of natural boundaries. This eigenvalue spectrum is shown 
schematically in figure [l] where the value of A con is shown by a horizontal dotted line, the eigenvalues 
Xn,L < A C on — by solid lines and the eigenvalues X Hi i 3= A con — by dashed lines. Let M(L) be the maximal 
value of n for which X„, L < A con and k n - M (L),L = [2{X n ,L ~ X con ) ID] 1 2 for n > M{L) and y/ c ° n (x) = 
Vn.ax) for n > M{L). Hence, we have 



, Y{x) M ^ , , t 

pax.t) = i — y e An ' Lt ys n ,ax )wn,ax) 

w Y{Xq) £ 



+ \l VTT e " Acon ' E e4Dfc ™' L Vfc " ^o)<" (x) . (2.27) 



m=l 

The solution with natural boundaries is the limit case L^oo 



p(x, t) = lim pax, t) (2.28) 

L— CO 
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V s (x) 




L72 X 



Figure 1. A schematic view of the eigenvalue spectrum for the problem with two absorbing boundaries 
at x = +LI2. The Schrodinger potential V${x) together with the boundaries at x = +LI2 is indicated by a 
solid curve and vertical lines. 



or 



Y(x) N ~ l 

p(X,t) = J— — Y. e nt Vn(Xo)Wn(x) 
I (Xq) „ =0 



Y { y\ 00 l 2 

li^e-Acon t Um £ e - 2 Dk mL t con (jCo) ^con {x) (2 2g) 

Y[x ) L -°°„T^i '"■ L rai 



where N = lim^oo M(L) is the number of bounded states in the case with natural boundaries. Since the 
eigenfunctions cannot be normalized at L — «■ 00, it is appropriate to write equation <2.29> for unnormal- 
ized eigenfunctions wi on (x), 



plx.t) = \ —— £ e Ant Vn(x )y/ n {x) 

Y(X ) n=Q 



F(x ) A k L " hL m ' L 



g-Hk,L) 

where the normalization constant JV is given by 

112 



= J dx|^ c fc on (x)| 2 (2.31) 



-LI2 

and the expression under infinite sum is divided and multiplied by Afcf, = fc m +i,L - k m> L- 

The infinite sum can be split into two parts: one with odd m and the other with even m. If the 
Schrodinger potential is symmetric, then one of these two parts contains only odd eigenfunctions i//^(x), 
whereas the other part has only even eigenfunctions i/z^x). In the limit L — «■ 00, these two sums can be 
represented by corresponding integrals, yielding 



Y{x) N ~ l 

p{x,t) = J — — - £ e" A »Vnto)VnW 
Y(Xq) n=0 



Y{x) 
Y{Xo) 



Y{xo) 



™' jdke-z D0t g;Hk)f e k (xa)tl(x) 


00 

Yix) r 1 2 

■e-^ 1 ldke-2 Dk tg-HVflixotylix), (2.32) 
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where 



ge(fc) = Jim 

L— »oo 



go(fc) = lim 

L— »oo 



L/2 



2Afc L J dx|^(x)| 2 



-L/2 
112 



2Ak L J dx\y/° k (x)\ 

-L/2 



This representation is useful if the eigenvalues and eigenfunctions are known. 



(2.33) 



(2.34) 



3. The analytical solution of FPE with constant force 

Let us consider a constant force term. In this case, the Fokker-Planck equation H2.lt reads 

dp[x,t) dp(x,t) Dd 2 p[x,t) 
= drift : + ' 



dt 



dx 



dx 2 



(3.1) 



This is a drift-diffusion problem for the potential V {x) = - v^stX normalized to V[x = 0) = 0. No station- 
ary solution exists for this problem, because the normalization constant jV in equation 12. lit diverges in 
this case. Nevertheless, the transformation (2.12) p{x, t) = Y{x) ll2 q{x, t) with 



Y{x) = exp 



2 V{x) 



, fdrift 

= exp|-^-x 



D 2 

can be used here to obtain an equation of Schrbdinger type 42.13) with constant Schrbdinger potential 

1 



(3.2) 



2D 



y drift ■ 



The stationary Schrbdinger-type equation corresponding to 12.211 reads 



d 2 y„(x) 
dx 2 



^drift 2 

~b~ 2 ~ D n 



fnix) = 0. 



Let us now add two absorbing boundaries located at x = +L/2, where y/{x = +LI2) = 0. 

Only in the case of real k n = \2X n ID - v^^/D 2 ] 1 2 > equation i3A\ has non-trivial solutions 

i//, 3 (x) = A cos(fc„x) + _Bsin(fc„x) , 

which satisfy the boundary conditions. These solutions are 

yicos(fc„ iL x) 



Vn,L(x) = - 



where n = 0,1,2,... and 



sin(fc„, L x) 



k n ,L = -(n+1). 



if n is even, 
if n is odd, 



According to l[3J>)-j377), we have from (233i and ( t2~34t 

ge(fc) = go(fc) = n . 

Taking into account that 

' D 



A con = lim min{A„, L } = lim mm\ —k nL + 

L— oo ■ 1 ' ' 



L— oo 



"drift 

2D 



"drift 
2D 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 
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holds, we obtain from equation 12.321 the expression 

f 1 

p tx, t) = exp — i> dri f t (x - xq ) 





\ v 2 1 
drift t 


exp 


2D 



i r i/2 

x— |dfce~2 Dfc *[cos(fcx)cos(fcxo) + sin(fcx)sin(fcxo)]. (3.10) 
n J 





Using the well known identities 

cos(fcx) cos(fcxo) + sin(fcx) sin(fcxo) = cos[fc(x- xo)] 

and 



J dk e" afc2 cosOSA;) = J ^e"^' 4 " , 



after simplification we obtain the well known result 



p{x, t) ■ 



>2Dt 



: exp 



4a 



(X-X - Adrift tY 



2Dt 



which describes a moving and broadening Gaussian profile. 



(3.11) 



(3.12) 



(3.13) 



4. Fokker-Planck dynamics with Poschl-Teller potential 

Here, as a particular example we consider the force 

f{x) = -fotanh(ax) (4.1) 
with some positive constants b and a. This corresponds to the diffusion problem in the potential 

b 

V(x) = — In (cosh ax) , (4.2) 
a 

normalized to V(x = 0) = 0. 

Figure[2]shows that this potential is actually a smoothed version of the V-shaped potential. The corre- 
sponding Schrodinger potential in this case is 



b 2 j b 2 ba\ 1 
Vs(x) = — -I — + 



2D [2D 2 ) cosh z (ax) 



(4.3) 




Figure 2. Graphical representation of equation 14.21 for b = 1 and several values of parameter a 
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Figure 3. Poschl-Teller potential (4.4) for Vq = 1 and several values of the parameter a. 



If we compare it [see equation (4.4) and figure [3) with the well known Poschl-Teller potential 

b 2 V 



VftM = VsM- 



2D cosh 2 (ax) 



(4.4) 



we see that equation (4.3) represents the shifted by b 2 /2D Poschl-Teller potential with Vq = b 2 l2D+bal2. 
As we can see from figure[3] the Poschl-Teller potential gives a mixed (discrete and continuous) eigenvalue 
spectrum. Therefore, equation <2.25> cannot be directly used to solve the FPE. We have to use <2.32t . 
The eigenvalue equation <2.16t for the potential d4.3t reads 



D d 2 y/ n {x) 
2 Ax 2 



b 2 ba 
+ 



2D [2D 2 ) cosh 2 (ax) 



1 



y/„{x) = -A, n y/ n {x) . 



(4.5) 



By introducing dimensionless variables x = ax, / = bIDa and \ n = 2X n IDa 2 - I 2 , we write <4.5> in a 
dimensionless form 



d 2 y/ n {x) 



dx 



v2 



1(1+ 1) 



1 



cosh x 



-y/ n (.x) = X n y/ n {x). 



(4.6) 



Analytical solutions for both bounded and unbounded eigenfunctions of equation d4.6t are known and 
can be found in Ha Hal . 



4.1. Bounded solutions for Poschl-Teller potential 



The equation <4.6( has N = maxjm e M | m < 1 + 1} bounded states n = 0, l,2,...,iV-l, where N is a set 
of all natural numbers N = {0, 1,2,...}. Here, we consider the eigenfunctions with \ n — as unbounded, 
because they cannot be normalized. 

The eigenvalues can be calculated from the following equation (lill 

A„ = -(/-n) 2 , for n<N; raeN. (4.7) 

Note that at least one bounded state with Xq- - I 2 always exists for Z > 0, which corresponds to Ao = 0. 
The bounded eigenfunctions are known 11211 

r [ jY e (n) Ff — i n, \ n — I; h; — sinh 2 x) if n is even, 

w„(x) = cosh"'(x)x^ v 2 ' 2 ' 2 1 , -o , , (4.8) 

l^ («)sinh(x)F(i-f,f + i-/;|;-sinh 2 x) if « is odd, 

where F denotes a hypergeometric function, which can be represented by Gaussian hypergeometric se- 
ries 

F(a,p;y;<) = > . (4.9) 

T{a)T{p)£ a T{j+k) nl 
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The normalization constants are 

n 1/2 



JY {n) = 



2(1 -n) 



(l-\n)[n + l) B (i,r-i/i)B(i 1+in) 
2(1- n) 1 



1/2 



(4.10) 



(4.11) 



I- \{n+ 1) B(|,/ - \{n+\))B(\,\{n+\)) 
where B(a, b) is the beta function B{a,b) = T{a)T{b)IT(a+ b). 

4.2. Unbounded solutions for Poschl-Teller potential 

The unbounded solutions have a continuous eigenvalue spectrum with =S A < oo. Thus, we can in- 
troduce k = A 1/2 (with k = kla). The Poschl-Teller potential is symmetric. Therefore, the eigenfunctions 
are the even and odd functions known from 1 13] 

<fr u m = A- y e kI (x)+B- (x) , (4.12) 

.(x) = (coshx)' +1 F(r,s;i;-sinh 2 x), (4.13) 

y/^jix) = (coshx)' +1 sinh(x)F(r+ -,s+ i;^;-sinh 2 Jcj, (4.14) 
where A and B are constants, and 

r = i(/ + l + ijfc), s=i(/ + l-ifc) . (4.15) 

Since these are unbounded solutions, eigenfunctions cannot be normalized within x e (-oo; +00). 

As we see, the eigenfunctions are rather complicated in general case. The expressions become es- 
sentially simpler for integer values of I. Therefore, without loosing the general idea, we will show the 
solutions of the Fokker-Planck equation for I = 1 and 1-2. 

4.3. The solution of FPE for Poschl-Teller potential with parameter 1 = 1 

For I - I (which implies b = aD) we have only one bounded state with the eigenvalue Ao = - 1 and 
the eigenfunction [equation H4.8) for n = 0] 

y/ m = —-^ . (4.16) 

V2cosh(jc) 

The unbounded eigenfunctions 14.131 and 14.14> are 

(jc) = cos(fcx) - -I tanh(x) sin(fcx) , (4.17) 
k k 

xfr° (x) = sin(fcx) + -I tanh(x) cos(fcx) . (4.18) 
* k 

As proposed in section[2] we add two absorbing boundaries located at x = +LI2. Due to these bound- 
ary conditions, we have only discrete values of k. Let us denote them by k- L m for even functions and by 
i<l m for odd functions. The values of k~ L and m , obtained from the boundary conditions, are positive 
solutions of the transcendent equations 

k Lm = tanh(I/2) tan(fc £ m L/2), (4.19) 
K L m tan(x L m ) = - tanh(L/2), (4.20) 
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where m = 1, 2,3, . . . denotes the m-th smallest positive solution. The equations for normalized eigenfunc- 
tions now read as 



K L,m 



cos{ki m x) - tanh(x) sin(fc^ m x) 

^L,m 



sin(K£ m x) + tanh(x) cos()C£ m x] 

K i,m 



where normalization constants for odd and even eigenfunctions are 

(fc 2 + l)[fcl-sin(fcl)] 



jr e {k,L) = - 



2fc 3 

(P + l) [kL+ sin(ifcl)] 
2p 



(4.21) 
(4.22) 

(4.23) 
(4.24) 



In the limit case L ^ oo, equations <4.19Hl4~2"0l for the allowed k values, as well as equations H4.231 - 
( 14.24) for the normalization constants simplify to 



V L— oo,m 



2mn 

L ' 
{2m -l)n 



In 

^k t =—, 



jV e [k,L—> oo) = JV (fc, L — >■ oo) 



and we also have 



ge(fc) = ^i^ oa -J/' e (k,L^oo) = n- 



L 
In 

T 

Lk 2 + l 
2 k 2 

k 2 + l 



k 2 

K 2 + l 



g [it) = Ak^^ • J/ (k, L—>oo) = ;r— 



(4.25) 
(4.26) 
(4.27) 

(4.28) 
(4.29) 



Inserting these relations as well as A con = l 2 a 2 D/2 (following from A con = 2X coa IDa 2 - 1 2 = 0) into (232) , 
we finally obtain the time-dependent solution of the Fokker-Planck equation 



p{x, t) = 



1 



2 cosh 2 (ax) 

| cosh(axo) c _i Da i t 
7rcosh(ax) 



oo 


oo 



dke \ Dal ~ k2t Ji — w|(ax)w|(axo) 
k 2 + 1 k k 



+ cosh^ e4Da2t r d - e4l)a a Ft J^ 

7rcosh(ax) J k 2 + l k k 

o 



(4.30) 



If the initial condition is given by xq = 0, then y/? (0) = and ifr'Z (0) = 1 hold, which allows us to obtain a 
simpler expression 

p{x, t) = ; : e - 



Da £ t 



2 cosh 2 (ax) 7rcosh(ax) 



o 



±Da 2 Pt k 



k 2 + l 



cos(fcax) - — tanh(ax) sin(A;ax) 
k 



(4.31) 



The solution for parameters b = 2, D = 2 and a = 1, corresponding to I = 1, with the initial loca- 
tion of the delta-peak at xq = 5 is shown in figure [4] for different time moments t . As we can see, the 
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Figure 4. The probability distribution at different time moments t, calculated for the parameters b - 2, 
D-2 and a = 1 (Z = 1) starting at xq - 5. 



probability distribution moves to the left. It broadens at the beginning. For larger times, it becomes nar- 
rower again and converges to the stationary solution p s tM = lim^oo p{x, t) = [cosh 2 (ax)] -1 = y/o(x) 2 
[see equations 1 14.30) and d4.16H , which is a symmetric distribution around x — 0. The stationary solution 
is practically reached at t = 10. This behavior is expected from the drift-diffusion dynamics. 

For small times t — «■ 0, we have a delta-peak located at x - xq in accordance with the given initial 
condition (2.3) . For comparison, the "general solution" of 1 10] does not satisfy this initial condition due 
to a wrong construction, where the contribution of bounded states is simply summed up with a Gaus- 
sian probability density profile (calculated with an error). The latter corresponds to unbounded states 
for zero Schrodinger potential at L — ► oo, as it is evident from 1 13.13) and 1 13.3) at i^nf* - 0- Therefore, the 
result appears to be correct only at f — ► oo when the Gaussian part vanishes. It is clear that the whole 
set of eigenfunctions should be calculated self-consistently for the given potential to obtain a correct and 
meaningful result, since only in this case the completeness relation 12.20) holds and all different eigen- 
functions are orthogonal. Thus, the basic error of 1 10] is that some of the eigenfunctions are calculated 
for zero Schrodinger potential in 1 10], whereas all of them should be calculated for the true Schrodinger 
potential. 



4.4. The solution of FPE for Poschl-Teller potential with parameter 1 = 2 

For 1-2 (which implies b = 2aD) we have two bounded states with eigenvalues Ao = 
The corresponding eigenfunctions are 



-4 and A 



ifo(x) = 
Vi (x) = 



V3 



2cosh 2 (Jc) ' 
3 sinh(x) 



2 cosh (x) 



The unbounded eigenfunctions are 



3tanh 2 (x) I cos(fcjc) - 3fctanh(x) sin(fcx) , 



1 + fc z - 3tanh 2 (x) I sin(fcx) + 3fctanh(x) cos(fcx) . 



i = -l. 

(4.32) 
(4.33) 



(4.34) 
(4.3S) 



By adding again two absorbing boundaries at x — +L/2, we have discrete values of k, i. e., ki m for 
even functions and i<i m for odd functions. In the limit L — «■ oo, we again obtain the classical infinite- 
square-well relations for eigenstates: 



1r~ 

^L— oo,m 



L^co,m 



(2m- 1)71 



2m7T 



(4.36) 
(4.37) 
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Figure 5. The probability distribution at different time moments t, calculated for the parameters b = 2, 
D = 4 and a = 1 (Z = 2) starting at xq = 5. 



The normalization constants in this case are 

^(fc.i — oo) = jY [k,L^ oo) = ^(P + 4)(fc 2 + l). 

By applying the same steps as in the case of I = 1, we obtain the solution 

3 6 sinniaxi sinmaxn) 3 n „2, 



p{x, f) = 



3 sinh(ax) sinh(axn) _j, 



4 cosh 4 (ax) 2 cosh 4 (ax) 
cosh 2 (ax ) -2Da 2 t 



n cosh (ax) 

,2, 



| cosh (ax ) c _ 2Da 2 t | d ^ c i/; c 
7i cosh 2 (ax) 





DC 

7 



dfce 2 



12 y \ [ax) v\ {axo) 

k 2 + 5k z + 4 " K 



FTiFTI <.(ax)^(ax ). 



(4.38) 



(4.39) 



The solution for parameters b — 4, D — 2 and a — 1, corresponding to / = 2, with the initial condition 
given by xo = 5 is shown in figure[5]for different time moments t . The evolution of the probability distri- 
bution is very similar to that one shown in figure |4]for 1 = 1, with the only essential difference that the 
dynamics is faster and the distribution is somewhat narrower due to a deeper potential well. 



5. Conclusions 

Using the analogy of the Fokker-Planck equation with the Schrodinger equation, it has been shown 
how the time-dependent solution can be constructed in the case of mixed eigenvalue spectrum with 
free and bounded states. The method is based on the idea of introducing two absorbing boundaries at 
x = +L/2, considering the limit L — ► oo afterwards. Although this idea is similar to the one proposed ear- 
lier in IIOH . it is obvious that the problem is quite non-trivial, so that the oversimplified (i.e., erroneous) 
approach of (3 cannot be used — see discussion in the end of section l4~3l Analytical solutions have been 
found and analyzed in two examples of the Schrodinger potential being constant (constant force) and a 
shifted Poschl-Teller potential. For the latter potential, the analytical solutions have been compared with 
the results of the Crank-Nicolson numerical integration method, and the agreement within an error of 
10~ 7 has been found. The time evolution of the calculated probability distribution in these examples is 
consistent with the usual drift-diffusion dynamics. 
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fll< p03B'fl3aTI/1 piBHflHHfl ^OKKepa-ll/iaHKa, BI/IKOpi/ICTOBytOHI/l 

cnexTp 3MiujaHi/ix B/iacHi/ix 3HaneHb? 

M. Epi/iLjcPl fl. Kayny3cP, P. MaHK^ 

IHCTHTVT (|)i3HKW, YHiBepCWTeT M. POCTOK, D-18051 POCTOK, HiMeHHUHa 

Ihctwtvt MaTeMaTHKH i KOMn'KyrepHwx HayK, /laTBMCbKHfi yHiBepcuTeT, LV-1459 Pwra, /laTBifl 

AHa/iorifl piBHflHHfl cpoKKepa-ri/iaHKa (FPE) 3 piBHAHHAM LUpeflUHrepa fl03BO/ifle BWKOpwcTa™ MeTOfl KBaHTOBOi 
MexaHiKH fl/ifl 3HaxoflxeHHfl aHa/ii™4Horo po3B'fl3Ky FPE fl/ifl hh3kh Buna^KiB. npoTe, nonepeflHi floaiiflxeHHfl 
o6Me>KyBa^MCfl noTeHuja/iOM LUpeflUHrepa 3 flucKpeTHUM cneKTpoM B/iacHnx 3Ha4eHb. TyT mm noKaxeMO, ak 

L\ei/\ niflXifl MOXHa T3KOX 3aCTOCyBaTH flO CneKTpy 3MiLUaHWX B/iaCHHX 3Ha4eHb 3i 3B'fl3aHWMH i Bi/lbHUMH CTa- 

HaMH. Mv\ po3B'fl3yeMO FPE 3 rpaHni^AMM, U40 3HaxoflflTbCfl npn x = +LI2 i 6epeMO rpaHnuro L — > 00, po3Mfl- 
flaKDHH npwiaiaflu 3 nocTMHWM noTeHuja/iOM LUpeflWHrepa i noTaHuja/iOM new/ifl-Te/i/iepa. Cnpom.eHwi'i niflxifl 
paHiiue 3anponoHyBa/iw M.T. Apayxo Ta E. flpiro Oi/ibfio. fleia/ibHe ,qoaii,qxeHHfl flBox npwuiafliB noKa3ye, U40 
KopeKTHMM po3B'fl30K, OTpi/iMaHuft b u,\v\ CTaTTi, y3roflxy£TbCfl 3 OMiKyBaHOK) flWHaiviiKOK) <t>OKKepa-n/iaHKa. 

KnioMOBi csioBa: pibhuhhh 0OKKepa-ll/iaHKa, pibhrhhr LUpefltiHrepa, noTeHu,iaji neuina-Tennepa 
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